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We study the effect of electric charge in compact stars assuming that the charge distribution is proportional to 
the mass density. The pressure and the density of the matter inside the stars are large, and the gravitational field 
is intense. This indicates that electric charge and a strong electric field can also be present. The relativistic hy- 
drostatic equilibrium equation, i.e., the Tolman-Oppenheimer-Volkoff equation, is modified in order to include 
electric charge. We perform a detailed numerical study of the effect of electric charge using a polytropic equa- 
tion of state. We conclude that in order to see any appreciable effect on the phenomenology of the compact stars, 
the electric fields have to be huge (~ 10 21 V/m), which implies that the total charge is Q ~ 10 20 Coulomb. 
From the local effect of the forces experienced on a single charged particle, it is expected that each individual 
charged particle is quickly ejected from the star. This in turn produces a huge force imbalance, and the gravita- 
tional force overwhelms the repulsive Coulomb and fluid pressure forces. The star can then collapse to form a 
charged black hole before all the charge leaves the system. 

PACS numbers: 04.20.Cv - 04.40.Dg - 04.40.Nr - 95.30.Sf - 97.10.Q - 97.10.Cv - 97.60.Jd - 97.60.Sm 



I. INTRODUCTION 

The study of the effect of electric charge and electric field 
in a gravitationally bound system has been done previously 
by some authors. Rosseland, in 1924 1 1] (see also Eddington 
01), studied the possibility that a self gravitating star on Ed- 
dington's theory, where the star is modeled by a ball of hot 
ionized gas, did contain a net charge. In such a system the 
electrons (lighter particles) tend to rise to the top because of 
the difference in the partial pressure of electrons compared to 
that of ions (heavier particles). The motion of electrons to the 
top and further escape from the star is stopped by the elec- 
tric field created by the charge separation. The equilibrium 
is attained after some amount of electrons escape leaving be- 
hind an electrified star whose net positive charge is of about 
100 Coulomb per solar mass, and building an interstellar gas 
with a net negative charge. As shown by Bally and Harri- 
son 0|, this result applies to any bound system whose size 
is smaller than the Debye length of the surrounding media. In 
fact, one should expect a star like the sun to hold some amount 
of net charge due to the much more frequent escape of elec- 
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trons than that of protons. Moreover, one should also expect 
that the escape would stop when the electrostatic energy of 
an electron e$ is of the order of its thermal energy kT. This 
gives for a ball of hot matter with the sun radius, a net charge 
Q ~ 6.7 x 10~ 6 T (in Coulomb). Hence, the escape effect 
cannot lead to a large net electric charge, and the conclusion 
is that a star formed by an initially neutral gas cannot acquire 
a net electric charge larger than about 100C per solar mass. 

In the case of cold stars, the gravitational pull is balanced by 
the degeneracy pressure of the particles. One can ask whether 
a net charge can have any effect on the structure of the sys- 
tem. For a star of mass M and charge Q, the electrostatic 
energy of a particle at radius r, eQ/r, is balanced by its grav- 
itational energy, mM/r, where e and m are the charge and 
mass of the particle, say a proton. For a spherical ball this 
gives a charge of approximately 100C per solar masses (see 
Glendenning|4|). This is clearly so for Newtonian stars. For 
compact stars, the high density and relativistic effects must be 
taken into account in order to reproduce with precision, the 
phenomenologies such as the mass and the radius. These ef- 
fects also reflect, in principle, in the allowed net charge of 
a compact star, the star can take some more charge to be in 
equilibrium (see e.g. Bekenstein |J|). 

That the observed stars, like the Sun or a neutron star, in 
equilibrium cannot support a great amount of charge comes 
from the fact that the particles that compose a star have a huge 
charge to mass ratio, as is the case for a proton or an electron. 
For highly relativistic stars, full general relativity is necessary, 
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and the situation might be different. One can have highly com- 
pact stars, whose radius is on the verge of forming an event 
horizon, such that the huge gravitational pull can be balanced 
by huge amounts of net charge. This type of configurations 
were raised by Bekenstein |5] and further developed in the 
studies by Zhang et al. @, de Felice and Yu fzl, de Felice et 
al. |8], Yu and Liu |9], Anninos and Rothman 1 10] and others. 
For instance, Zhang et al. 1 6] found that the structure of a neu- 
tron star, for a degenerate relativistic fermi gas, is significantly 
affected by the electric charge just when the charge density is 
close to the mass density (in geometric units). In the inves- 
tigations by de Felice et al., and by Anninos and Rothman, 
they assumed that the charge distribution followed particular 
functions of the radial coordinate. 

In passing, we can mention that if the charge to mass ra- 
tio of the particles that make the star is low, say one or of 
order one, then a star can contain a huge amount of charge. 
For instance, one could think of dust particles containing a 
large quantity of neutral particles, say 10 18 neutrons, and one 
proton. This would give a charge to mass ratio of the order 
one. The dust particles could then cluster around each other 
and form a star. Theoretically these stars have been consid- 
ered. The Newtonian theory of gravitation and electrostatics 
admit equilibrium configurations of charged fluids where the 
charge density p c h can be as large as the mass density p, in 
appropriate units. For instance, a static continuous distribu- 
tion of charged dust matter with zero pressure will be in equi- 
librium everywhere if p c h = ±p (in geometric units). This 
gives for a ball of dust with the mass of the sun a net charge 
of about 1.7 x 10 20 Coulomb. The general relativistic ana- 
log for charged dust stars was discovered by Majumdar 1 11] 
and by Papapetrou 1 12], and further discussed by Bonnor 11311 
and several other authors (see 1 14] for a review). We will not 
pursue this line here. 

In this paper we shall take the ideas initiated in |5] and con- 
tinue in oBOiElQlJL an d focus attention on the effect of 
charge in cold compact stars, made of neutrons, protons and 
electrons, say. We assume that the charge density goes with 
mass density and write p c h = ap, a being a positive (neg- 
ative) constant for a positive (negative) charge density. We 
further assume that the net charge in the system is in the form 
of trapped charged particles carrying positive electric charge. 
Later we will show that the negative charge will also give the 
same results, thus the sign of the charge has no effects. From 
the Einstein-Maxwell equations, we write a modified Tolman- 
Oppenheimer-Volkoff (TOV) equation (see 01) where the en- 
ergy density which appears from the electrostatic field will 
add up to the total energy density of the system, which in turn 
will help in the gaining of the total mass of the system. The 
fundamental difference to the standard relativistic TOV equa- 
tion from the uncharged case is the presence of the Maxwell 
stress tensor, which in the modified TOV equation contributes 
to the energy density, and to the pressure, besides the pres- 
ence of a potential gradient term from the coulombian term. 
We solve the modified TOV equation for polytropic equations 
of state assuming that the charge density goes with the matter 
density and discuss the results. We will show that a previous 
estimate by Bekenstein is correct jjt : in order to see any ap- 



preciable effect on the phenomenology of the neutron stars, 
the charge and the electrical fields have to be huge, the large 
electric field being able to start producing pair creation of par- 
ticles. 

Furthermore, we find solutions with huge amounts of 
charge. How this extra charge is formed in the star is not the 
concern here. A mechanism to generate charge asymmetry 
for charged black holes has been suggested recently 11511 and 
the same may be applied for compact stars too. We are not 
claiming that compact stars always have so large a charge and 
strong electric field. Indeed, one has to raise the question what 
happens to a single charged particle inside the system. It will 
face a huge electrostatic force and will soon escape the sys- 
tem making the star unstable. During that process, the gravi- 
tational force, which was previously balanced to some extent 
by the repulsive Coulombian force, will overpower the inside 
pressure and the system will collapse to form a black hole. 
However, not all the charge has time to escape from the sys- 
tem and hence the charge will remain trapped inside the event 
horizon to form a charged black hole. This scenario could 
be applicable for newly born charged compact stars and their 
rapid transformation to charged black holes. 

This paper is organized in the following sections. In section 
im we give the general relativistic formulation for the inclu- 
sion of charge in the TOV equation and present the numerical 
procedure to solve them for a given radial charge distribution 
inside the star. Section presents the application of these 
equations to polytropic stars wherein different aspects of the 
effect of charge are described, such as mass-radius-charge di- 
agrams, the inside electric fields and the metric. In section 
II VI we discuss the stability of these stars and the possibility 
of formation of charged black holes. Finally we draw our con- 
clusions in sectionM 

II. GENERAL RELATIVISTIC FORMULATION 

We take the metric for our static spherical star as 

ds 2 = e v c 2 dt 2 - e x dr 2 - r 2 (d6 2 + sin 2 9d<j, 2 ). (1) 

The stress tensor will include the terms from the 
Maxwell's equation and the complete form of the Einstein- 
Maxwell stress tensor will be : 

Tj} = (P + eK^-P^ + i- (V a F Q „ - -^F aP F a ^j 

(2) 

where P is the pressure, e is the energy density (=pc 2 ) and 
is 4- velocity vector. For the time component, one easily sees 
that u* = e~ v ' 2 /c and hence u*it t = 1. Consequently, the 
other components (radial and spherical) of the four vector are 
absent. 

Now, the electromagnetic field is taken from the Maxwell's 
field equations and hence they will follow the relation 

[V~9F^] v = 4Kj"V=g (3) 

where j M is the four-current density. Since the present choice 
of the electromagnetic field is only due to charge, we have 
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only F 01 = — F 10 , and the other terms are absent. In gen- 
eral, we can derive the electromagnetic field tensor f" from 
the four-potential A^. So, for non vanishing field tensor, the 
surviving potential is Aq = <fi. We also considered that the 
potential has a spherical symmetry, i.e., tj> = </>(r). 

The non vanishing term in Eq. (0 is when v = r. This 
gives the electric field as 

F» a F au - hnF a0 F a ^ - ~ U ~ 
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where, 



U{ r ) = 1 f ^r 2 fe^' 2 dr 
r Jo 



(4) 



is the electric field. Let us consider a distribution of charge in 
the matter as some charge density annoted as p^. We modify 
the above expression (Eq. and rewrite it as 



1 f r 

U(r) = — 4Trr 2 p ch e x/2 dr. 
r Jo 

he total charge of the 

/ 4Trr 2 p ch e x/2 dr 
Jo 



(5) 



This leads to defining the total charge of the system as 

Q= I 4nr z p ch e A/z dr (6) 



where R is the radius of the star. 

With the metric Q one can easily get the Einstein's field 
equations from the relation 
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The first of the Einstein's equations is used to determine the 
metric e A . The mass of the star is now due to the total con- 
tribution of the energy density of the matter and the electric 
energy (^) density. The mass takes the new form as 
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and the metric coefficient is given by 

2GM„(r) 
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However, this mass is the one measured in the star's frame. 
For an observer at infinity, the mass OlQfJ ls given by 

U 2 



Moo = 



4irr 



inr 2 



R 

M tot (R) - 



8ttc 2 
U 2 
8ttc 2 
U 2 



cr 

Q(R) 

2R 



8ttc 2 / 



dr 
dr + 
dr 



(11) 



where R is the radius of the star. In our plots and rest of the 
text, we will refer to the mass M as the M to t- 

From the conservation of stress tensor (T^ = 0) one gets 
the form of the hydrostatic equation. To this end, we obtain 
the modified TOV as 
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The first term on the r.h.s. comes form the gravitational force 
with an effective pressure and density which we will discuss 
later, and the second term from the Coulomb force that de- 
pends on the matter by the metric coefficient. 

Numerical solutions in fortran require the integral forms of 
Eqs. (0, ® & d 1 0I> to be expressed in their differential forms. 
We can write the corresponding differential forms as 
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and 
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U 2 
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dr. (15) 



So, the final four equations needed to be solved are Eqs. \\2\ . 
dl3> . dl4> & (I15> . The boundary conditions for the solution 

are, at the centre where r = 0, U(r) — 0, e A ( r ' = 1, P{r) = 
P c and p{r) — p c , and at the surface where r = R, P(r) = 0. 
The inputs in the equations are the pressure P, the energy 
density e and the charge density p c ^. The metric coefficient 
A and the electric field U are interdependent. This gives us 
a set of four coupled differential equations which we solve 
simultaneously to get our results. We note here that the form 
of the equations does not change with the sign of the charge 
because the electric field appears in the mass term 14l > and 
pressure gradient term dl2l in squares and in the Coulomb 
part, the product p c \JJ is also invariant. 



III. CHARGED STARS IN A POLYTROPIC EQUATION OF 

STATE 

A. The Mass-Radius relation and other features 

We now examine the effects of charge in a polytropic equa- 
tion of state (EOS), which is a more general approach than 
considering any model dependent EOS. We make the charge 
go with the mass density (e) as 



Pch = f X 



(16) 



where e = pc 2 is in [MeV/fm 3 ]. With this assumption, the 
charge fraction / has a dimension ^^[frr^/MeV] 1 / 2 and 
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the charge density p c h is in [MeV/fm 3 ] 1 / 2 -!^^. This kind 
of 'mixed units' appear in our dimensions of / and /? c h be- 
cause one can see from Eqs. l|5} & d!2t that the electric field 
is proportional to the square root of the pressure (in units of 
MeV/fm 3 ) and the integration over the radius r is carried out 
in kilometers. In geometrical units, this can be written as 



p ch = a x p 



(17) 



where charge is expressed in units of mass and charge den- 
sity in units of mass density. This a is related to our charge 
fraction / as 



a = f x 



0.224536 



G 



f x 0.86924 x 10 3 



(18) 



Our choice of charge distribution is a reasonable assumption 
in the sense that large mass can hold a large amount of charge. 
The polytropic EOS is given by 



P 



Kp 



1+1/ra 



(19) 



where n is the polytropic index and is related to the exponent 
r as n = p^rj". In the relativistic regime, the allowed value of 
T is | to | . We have considered the case of T = | and the 
corresponding value of n is 1.5. Primarily, our units of matter 
density and pressure are in MeV/fm 3 . We choose a value of 
k as 0.05 [/m] 8 / 3 . Thus we have an equation of state which 
we analyze for different cases of charge fraction / and study 
the nature and behavior of the system. The choice of the EOS 
is however not very important in so far as the nature of the 
curves due to the effect of charge is concerned, and show only 
corresponding shift in the maximum mass based on the type 
of EOS used. 

It should be noted that the nuclear forces are only affected 
by electromagnetic forces (and so the EOS) when the number 
density of the charged particles are of the order of the baryonic 
number density, i.e., Z ~ A. This has been verified previously 
in many works on the nuclear structure. However, for the case 
of charged stars, the forces are compared between the elec- 
tric and the gravitational force. For our case, the Z/A ratio is 
~ 10~ 18 , which essentially means that these extra charged 
particles which produces huge electric field and affects the 
structure of the star, produces negligible effect on the nuclear 
matter and the EOS. The same argument holds for the chem- 
ical potential, which are controlled by chemical equilibrium 
and charge neutrality. The extra charged particles are just very 
few when compared to the total number of baryons of the star. 
Any effect brought in by these extra charged particles can be 
approximated (and compared) to the case of chemical equilib- 
rium and zero charge. So, it is justified to use an EOS which is 
calculated for neutral matter. We have indeed tested the effect 
of this amount on charge on the EOS on some sophisticated 
models of nuclear matter like the Non-linear Walecka model 
and verified our reasoning to be right. 

We plot the mass as a function of the central density p c in 
Fig. Q for different values of the charge fraction /. The stars 
which are on the higher density regime and have lower mass, 
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FIG. 1 : Central density against mass for different values of the factor 
/• 



are unstable because i^- < 0. Those falling in the lower 
density regime and have increasing mass are all allowed. The 
effect of the charge for / = 0.0001 on the structure of the star 
is not profound and is comparable with that of chargeless star. 
This value of / is however 'critical' in the sense that further 
increase in the value shows effect on the structure of the star. 
With the increase of the charge fraction form / = 0.0001 
to / = 0.0005 the structure changes by 20% increase in the 
value, from / = 0.0005 to / = 0.0008, the increase is 35% 
and from / = 0.0008 to / = 0.001, the change is almost 90%, 
thus showing that the change in the structure is non-linear with 
the change in the charge fraction as can be seen in Fig. Q. 
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FIG. 2: Mass-Radius relation for different values of the factor /. 

In Fig. (0 we plot the mass-radius relation. Due to the 
effect of the repulsive force, the charged stars have large ra- 
dius and larger mass as we should expect. Even if the radius 
is increasing with the mass, the M/R ratio is also increasing, 
but much slower. For the lower charge fractions, this increase 
in the radius is very small, but a look at the structure for the 
fraction / = 0.001 reveals that for a mass of 4.3 M©, the ra- 
dius goes as high as 35 km. Though the compactness of the 
stars are retained, they are now better to be called as 'compact 
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charged stars' rather than 'charged neutron stars'. The charge 
fraction in the limiting case of maximally allowed value goes 
up to / = 0.0011, for which the maximum mass stable star 
forms at a lower central density even smaller than the nuclear 
matter density. This extreme case is not shown in Fig. (0 be- 
cause the radius of the star and its mass is very high (68 km 
and 9.7 M Q respectively). These effects suppress the curves 
of the lower charge fractions due to scaling. For this star, the 
mass contribution from the electric energy density is 10% than 
that from the mass density. It can be checked by using rela- 
tion (1181 that this charge fraction / = 0.0011 corresponds to 
Pc\i = 0.95616 x p in geometrical units. In Table|I] we show 
the maximum mass for different charge fractions and their cor- 
responding radii, central densities and net charge content. 

TABLE I: The maximum allowed stable stellar configuration for dif- 
ferent charge fractions. 



/ 


Mtot 




R 


e c 


Q 




CM®) 


(M ) 


(km) 


(MeV/fm 3 ) 


(xlO 20 Coulomb) 


0.0001 


1.428 


1.43 


11.87 


1550.41 


0.259 


0.0005 


1.69 


1.765 


13.55 


1202.6 


1.517 


0.0008 


2.438 


2.728 


18.47 


652.87 


3.434 


0.001 


4.384 


5.248 


31.47 


226.55 


7.576 


0.0011 


9.76 


12.15 


68.7 


47.04 


18.314 



0.001 
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FIG. 3: Variation of the metric e with the radius in the maximum 
mass stars for different cases of charge distribution with varying /. 

In Fig. (0, we plot the metric coefficient e A as a func- 
tion of radius for the maximum mass stars for each of the 
charge fractions. A quick comparison shows that the na- 
ture of e A is the same for all the stars with different charges. 
There is a slight increase in its value for higher charge frac- 
tions thus showing the gain in the compactness of the star 
with charge. This can be verified from the values of e A at 
the surface for two cases of charge fractions / = 0.0001 
and / = 0.001, (^r) nnnm = 1 - tte = 0.3548, and 

■> ' V R 1 0.0001 1.55 

(tt)oooi = 1 _ IT = °- 4118 - This re-confirms that the 



compactness of the star increases despite the enormous in- 
crease in the radius with the increase of the charge fraction. 

It is interesting to point out that we have tested our model 
of charged stars for a very soft EOS which can go upto a very 
high density, where we saw that the spiraling behaviour of the 
mass-radius curve exists for high charges. This nature of the 
curve are very much supportive to the fact that these charged 
stars are stable to stellar oscillations. 



B. Charge and Electric field inside the star. 

In this subsection we discuss the effect of the charge in 
the inner profile of the star. In the Q x R diagram in Fig. 
@, we plot the total charge Q of the stars at the surface as 
function of their radius. From Eq. (J6J, it is clear that the 
charge will increase with the increasing charge fraction. This 
increase is however not the only one responsible for the high 
charges if we consider that equation. There are contributions 
also from the large radius of the higher charge configuration 
and the larger value of the metric coefficient e A . So, we see 
that the matter and charge are inter-related to each other and 
the effect of one depends intrinsically on the other and hence 
the relations are strongly coupled. 
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FIG. 4: The variation of the charge with radius for different /. 

The charge developed due to this range of charge fraction 
(/ = 0.0001 to / = 0.0011) is in the limiting case where the 
effect is noticed on the structure. This scale of the charge is 
easy to understand from the mass expression and the modified 
TOV (Eqs. (|9j and (I12» where in order to see any change in 
the stellar configuration, we should have 

U ~ V8wP < s/Sne. (20) 
If we consider that P ~ 10 MeV/fm 3 , then U ~ 

1 /2 

10 (^r) an d we show in appendix that with proper con- 
version, this gives U ~ 10 22 V/m. With R~10 km, the charge 
needs to be at the order of 10 20 Coulombs. 

In the Q x M diagram in Fig. (JSJl, we plot the mass of 
the stars against their surface charge. We have made the 
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FIG. 5: The variation of the charge with mass for different /. 



charge density proportional to the energy density and so it 
was expected that the charge, which is a volume integral of the 
charge density, will go in the same way as the mass, which is 
also a volume integral over the mass density. The slope of the 
curves comes from the different charge fractions. The nature 
of the curves in fact reflect that charge varies with mass (with 
the turning back of the curves all falling in the 'unstable zone' 
and is not taken into consideration). If we consider that the 
maximum charge allowed is estimated by the condition J20i 
for 4r- to be negative (Eq. H2\ ). we see that the curve for the 
maximum charge in Fig. (|5} has a slope of 1:1 (in a charge 
scale of 10 20 Coulomb). This scale can be easily understood 
if we write the charge as 



Q = VGMq— ~ 10 20 — Coulomb. (21) 

Mq Mq 

It is worth mentioning that this charge Q is the charge at the 
surface of the star where the pressure and also 4^ are already 
very small (ideally zero). So, at the surface, the Coulomb 
force is essentially balanced by the gravitational force and the 
relation of the charge and mass distribution we found is of the 
same order (i.e., Q ~ M or p c h — p) for the case of charged 
dust sphere discussed earlier by Papapetrou [12] and Bonnor 
iTnll . as we referred in the introduction. 

In Fig. 0, the electric field (U) is plotted as a function 
of the stellar radius for the maximum mass star for different 
charge fraction /. The radius is shown in a log-scale. From 
the expression (Eq. (0) of the electric field, it is clear that the 
profile of the field will depend on the charge fraction /, the 
metric coefficient e A and also the radius of the star. The value 
of the field increases for charge fraction up to / = 0.0008 
but then falls down. The decrease in the field for the higher 
charge fraction is attributed to the formation of the stable star 
at a very small density (Fig. Q) for that particular charge 
fraction /. 
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FIG. 6: Variation of U with the radius in the maximum mass stars for 
different cases of charge distribution with varying /. 

C. The modified EOS inside the stars 

^From Eqs. © and dl2b . the effective energy density and 
pressure of the system can be viewed as 

7/2 

P* = P (22) 

and 



In Figs. and (jSJl we show two different charge fractions 
/ = 0.0005 and / = 0.0008, the effective pressure (P*) and 
energy density (e*). The effective pressure drops down to a 
negative value, but even in this case, the first part of Eq. (II 2b 
preserves its overall negative value because (M to t + 47rr- 3 P*) 
is positive. Thus the overall sign of the pressure gradient 4p 
is still negative as long as the attractive gravitational term is 
larger than the repulsive Coulombian one. These figures show 
that the effective EOS becomes stiffer due to the inclusion of 
charge and consequently allowing more mass in the star. 

Also, the effective pressure directly reduces the value of 
the negative part of 4^, but the effective energy density in- 
creases the same through the A/ to t- This goes on until the ef- 
fective pressure becomes so much negative that it overcomes 
the value of the M to t and this limits the formation of star with 
higher charge fraction /. 

D. Balance between the Coulomb and gravitational forces 

In this section, we discuss the effects brought in the pres- 
sure gradient from the matter energy and the Coulomb energy. 
As mentioned previously, the total mass of the system M tot 
increases with increasing charge because the electric energy 
density 'adds on' to the mass energy density. This change 
in the mass is low for smaller charge fraction and going up to 
seven times the value of chargeless case for maximum allowed 
charge fraction / = 0.0011. This effect however does not 
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FIG. 7: The pressure versus energy density plot for chargeless case 
and for charge fraction / = 0.0005. 
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FIG. 9: The Pressure gradient with the positive Coulomb part and the 
gravitational negative part is shown here for two different values of 
the charge factor / (=0.0005 and 0.0008). The total, jZ, is also plot- 
ted for the same. The annotations for the component of pressure for 
two charges / = 0.0005 and 0.0008 coming from the matter part 
are dPg5 and dPg8 respectively, those from Coulomb part are dPc5 
and dPc8 respectively and the corresponding totals are dP5 and dP8. 

part and a matter part, both are coupled with each other in 
the sense that charge changes the mass and the mass changes 
the Coulomb part. The plots are for two values of the charge 
fraction / = 0.0005 and / = 0.0008. The positive part of ^ 
maintains its almost constant value because the charge frac- 
tion / is the controller of the same, and in our case, they differ 
by a very small percentage. In the negative part, the changes 
are drastic and are mainly brought by the effective pressure as 
we already discussed. 



FIG. 8: The pressure versus energy density plot for chargeless case 
and for charge fraction / = 0.0008. 



change the metric coefficient considerably when compared to 
the chargeless case (see Fig. Q). So, the ruling term in Eq. 
O is the factor (M tot + Anr 3 (P - ^)). With the increase 

of charge, the value of (P-^) decreases, and hence the gravi- 
tational negative part of Eq. H2\ decreases. The central pres- 
sure is very high. So, with the softening of the pressure gradi- 
ent, the system allows more radius for the star until it reaches 
the surface where the pressure (and the pressure gradient too) 

7/2 

goes to zero. We should stress that because cannot be too 
much larger than the pressure in order to maintain 4^ nega- 
tive as discussed before, so we have a limit on the charge. It 
is interesting that this limit comes from the relativistic effects 
of the gravitational force and not just only from the repulsive 
Coulombian part. 

This effect is illustrated in Fig. where we plot both the 
positive Coulomb part and the negative matter part of the pres- 
sure gradient. Note that, although we mention a Coulomb 



IV. STABILITY OF CHARGED STARS AND FORMATION 
OF CHARGED BLACK HOLES 

A. Stability considerations 

Here we will discuss the stability of these charged compact 
stars taking into account the forces acting only on the charged 
particles. For this, we need to compare our results with the 
stable configuration of an almost neutral neutron star. 

The basic argument to assume the charge neutrality of New- 
tonian stars, or quasi-Newtonian stars (such as neutron stars), 
is based on the fact that the total charge of the stars should 
lie below a certain limit where the Coulomb repulsive force 
overwhelms the gravitational attractive force at the surface of 
the star. This limiting value can be viewed as the Coulomb 
force acting on a proton at the surface and we have the limit- 
ing range as : 



(Ze)e GMm < G{Am)m 



R 2 



R 2 



R 2 



(24) 



where, Ze, R and M are the net charge of the star, its radius 
and mass respectively and m and e are the mass and charge 
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of a proton. In the above equation, the mass of the star is 
considered to be smaller than Am because of the gravitational 
binding of the system. In gravitational units, 

!<(?) 2 <10- (25, 

If we take into account that there are approximately A ~ 
10 57 baryons in a neutron star (M ~ M ), so Z < 10 21 , 
which gives 

Q ~ Ze ~ 100 Coulomb. (26) 

This is the limit on the net positive charge already discussed 
in the introduction. If the star had a net negative charge, then 
the electrons being the carrier of charge, the value of the limit 
would be reduced by the factor m e /m |4[. Using the solution 
for the TOV equation in the presence of electric fields, which 
we already solved in the previous sections, we showed that the 
balance of the forces, to make an element of the fluid at rest 
(hydrostatic equilibrium equation) allows a very large amount 
of charge. 

In fact we showed that the maximum allowed charge is ob- 
tained when p c h — VGp in natural units from the assumption 
that more mass can hold more charge (p c h oc p). Because 
of the large density found in the compact stars like neutron 
stars (p ~ 10 15 g/cm 3 ) we can expect a large charge if p c h is 
of the same order as p. Thus, Q(r) ~ VGM(r) and at the 
surface of a neutron star with M ~ M , the total charge is 
Q ~ %/GMq ~ 10 20 Coulomb and the intensity of the elec- 
tric field is of the order of 10 21 — 10 22 V/m. In the case of 
maximum charge, the relation Q(r) = \GM{r) at the sur- 
face shows that the ratio Zj A ~ 10 -18 as compared to 10 -36 
in Eq. J25> . This explains why we are having charges 10 18 
times larger than that in relation J26t . This Z is the net charge 
in the star and is the difference of the charged particles with 
opposite sign (Z ~ Z net ~ \Z + — Z-\). The number of 
charges of equal and opposite sign is not at all limited. 

However, not only a global balance of all the forces is 
enough to guarantee the stability of these stars when we con- 
sider the forces acting on each of the charged particles of that 
fluid. As the Z/A ratio is ~ 10~ 18 , the Coulomb force is 10 18 
times larger than the gravitational one felt by a charged par- 
ticle. This will make the charged particles leave the system, 
making the star unstable unless any other mechanism exists to 
bind together one charged particle with 10 18 neutral particles. 

Our analysis reveals that the maximum amount of charge 
which can be allowed in the system of charged compact stars 
is close to unity (precisely 0.95) in geometrical units and can 
be compared to the extremum case of charge in a Reissner- 
Nordstrom black hole, i.e., Q/y/GM = 1. The stability anal- 
ysis of the charged star from the viewpoint of the force acting 
on a single charged particle makes the star highly unstable, 
with the charge escaping from the star within a very short 
time. However, very little amount of charge will escape be- 
fore it collapses to a charged black hole. This can be seen as 
follows : with the escape of little amount of charge, the re- 
pulsive pressure from the Coulomb part diminishes, but the 
attractive pressure from the matter part is not affected very 



much. This is because the gravitational contribution of the es- 
caped charged particle (say proton) is practically zero as com- 
pared to that from the 10 18 neutral particles. This leads to 
a very critical scenario of disbalance of the global forces in- 
side the star, where the loss of a little amount of charge makes 
the gravitational attractive pressure overwhelm the residual re- 
pulsive Coulombian pressure and the whole system collapses 
further to form black holes. But all the charge has not yet es- 
caped from the system and so, the residual charge gets trapped 
in the black hole forming the charged black holes. This pic- 
ture is best suitable for the extreme charged case, but can be 
true for smaller charged fractions also. We saw in Fig. Q, 
that the maximum stable mass for the highest allowed charge 
fraction (/ = 0.0011) is as high as 9.7 M and such high 
mass charged compact stars are more favorable to collapse to 
Reissner-Nordstrom black holes. 

B. Discharge time 

Since the electric force on a charged particle is much greater 
than the gravitational force one can neglect the gravitational 
force, and find that, from Newtonian theory, the equation of 
motion for a charged particle is 

d 2 r _ Qq 
m W~ r 2 

where m and q are the mass and charge of the particle, and 
Q is the net charge of the star. Thus, the lifetime to dis- 
charge the star fdischarge is given in Newtonian theory by 

idischargoNcwt ^ \\J gm/c* *f ' where we have used the re- 
lation Q ~ \[GM. For R ~ 10 km, GM/c 2 - 3 km and 
m/q ~ 10~ 18 , we have, t ~ 10~ 13 s. However this time 
scale is much shorter than the time a particle takes to traverse 
the radius of the star. This means that the charged particle 
acquires quickly a velocity near the speed of light and it will 
travel throughout the star with this speed. Thus, the discharg- 
ing time is of the order of the crossing time, i.e., 

^discharge — 

For a star with R ~ 10 km tdischarge — 10~ 5 s. 

We have not taken into account that the particle collides 
with other particles with a mean free path of the order of one 
fermi. However, the trajectory of the charged particle is not 
going to be a random walk trajectory, since the electric force 
is radial. Thus, this collision process might increase a bit the 
discharging time, but not a lot. A detailed evaluation is how- 
ever a very complicated process and is beyond the scope of 
the present paper. 

V. CONCLUSIONS 

Our analysis shows that the amount of charge contained in 
a dense system like a compact star can be very high and sev- 
eral orders of magnitude larger than those calculated by clas- 
sical balance of forces at the surface of a star. This amount 
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of charge mainly comes from the high density of the sys- 
tem since p c h oc p. We showed that the charge can be as 
high as 10 20 Coulomb to bring in any change in the mass- 
radius relation of the star, yet remaining stable as long as one 
considers only the hydrostatic equilibrium and the global bal- 
ance of forces and not considering the individual forces on 
the charged particles. In our study, we used a polytropic EOS 
for our compact star with a choice of parameters such that the 
system is close to the realistic neutron stars. We showed that 
in the critical limit of the charge contained in the system, the 
maximum mass stable star forms in a lower density regime, 
however compactness keeps on increasing. We have also stud- 
ied the change in the pressure gradient due to the effect of the 
charge contained in the system. It was expected from the clas- 
sical picture of forces involved, that the repulsive force of the 
charged particles will add up to the internal pressure of the 
system and the entire repulsive force will be balanced by the 
gravitational force of the system. However, we can see from 
Fig. (|9) that the contributions from the charged particles are 
helping to soften the gravity part of the pressure gradient, thus 
allowing more matter to stay in the system. The second term 
in the right hand side of Eq. (I12t is positive always and does 
not depend on the nature of the charge whether it is positive or 
negative. The pressure gradient must be a negative term and 
hence 4? is softened by the effect of the presence of charge. 
The net effect is the gravitational force which tries to collapse 
the system is held further away by the Coulomb force, which, 
in the absence of the gravity would have exploded the system. 

This imbalance of the pressure actually will come to play 
when some of the charge will leave the star due to the Coulom- 
bian repulsion acting on a single charged particle. As we 
discussed in section IIVI the attractive gravitational pressure 
will then become more than the repulsive total kinetic and 
Coulombian pressure, as a result of which the star may col- 
lapse to a charged black hole. 

We have studied the electric field inside the maximum mass 
star allowed by a certain charge fraction. We found that the 
field attains a maximum value for certain amount of charge 
fraction and then decreases for higher charge fractions (Fig. 
0). This is also interesting because normally it was expected 
that the increase in the charge fraction would increase the field 
also. But as our charge distribution varies directly with the 
matter density, so the formation of the maximum mass star 
in the lower density regime for high charge fractions, reduces 
the electric field. The mass of the stars (Eq. l|9}) in these crit- 
ical field limit is very high as compared to the chargeless case 
and have a contribution from the electrostatic energy density, 
but this is small compared to the contribution of the matter 
density. 

As pointed by Bekenstein |5], with this amount of huge 
charge (as we have in our system), the electric field produced 
will be too high and give rise to pair production. This effect 
results in self-diminishing the electric field and thus the sys- 
tem will destabilize 11911 . However, we can say that the critical 
field limit calculation has been done in vacuum, and it is not at 
all clear how the field will behave in a dense system. Recent 
observations have revealed that there are magnetars, stars with 
strong rotating fields, which have a magnetic field as high as 
10 15 Gauss. The critical magnetic field limit for pair creation 



in vacuum is 10 13 Gauss. So, if such high fields really exist 
in a highly dense magnetar and these stars are stable, then the 
critical field limit needs to be modified for high density mat- 
ter. Putting aside this debate of the critical field limit, we have 
checked only the behaviour of the system in the presence of 
a very high charge. We have found that the global balance of 
forces between the matter part and the electrostatic part can 
support a huge amount of charge. In addition, we have argued 
that the solution is unstable, and a decrease in the electric field 
can create an enormous imbalance, resulting in the collapse of 
the charged star to a charged black hole. Finally, these charged 
stars should be very short lived, they would exist within the 
very short period between the supernova explosion and the 
formation of the charged black holes. 
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APPENDIX: CONVERSION TO REAL UNITS 

It is important to mention here the basic units used in our 
approach and their conversions from these units to the real 
units of the charge and fields like Coulomb and Volt/meter 
respectively. We used the charge density (p c h) proportional 
to the energy density of the system (p), which in turn is in 
MeV/fm 3 , with a factor / which has dimensional units. From 
the fine structure constant a 



he 



h=, we get a relation 



137 



1 statCoulomb(l e.s.u.) = 2.51 x 10 9 [MeV fm] 1 / 2 . Using 
the conversion factor 1 Coulomb = 2.998 x 10 9 c.s.u. for the 
charge we arrive at 



1 Coulomb ~ 0.75 x 10 19 [MeVfm] 



,1/2 



(A.l) 



u 2 

Also, the units of have to be the same as the units of 
pressure. Initially, we have that the units of pressure are 
[MeV/fm 3 ]. This indicates that U wiU have [MeV/fm 3 ] 1 / 2 
units. Working out the relations, we find 



MeV 



Also, 



fm 3 



MeV 



1/2 



= 1.2 x 10 21 V/m 



c 2 fm 3 



1.78 x 10 12 g/cm 3 



(A.2) 
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